Thermal collapse of spin-polarization in half- metallic ferromagnets 
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The temperature dependence of the magnetization and spin-polarization at the Fermi level is 
investigated for half-metallic ferromagnets. We reveal a new mechanism, where the hybridization of 
states forming the half-metallic gap depends on thermal spin fluctuations and the polarization can 
drop abruptly at temperatures much lower than the Curie point. We verify this for NiMnSb by ab 
initio calculations. The thermal properties are studied by mapping ab initio results to an extended 
Heisenberg model which includes longitudinal fluctuations and is solved by a Monte Carlo method. 
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Half-metallic ferromagnets (HMFs) are ferromagnetic 
metallic compounds showing, in the ideal case and at 
zero temperature, a spin polarization of P = 100% at 
the Fermi level Ep. This means that the spin-resolved 
density of states (DOS) shows a metallic character only 
for one spin direction (usually majority-spin), with en- 
ergy bands crossing Ep; contrary to this, the other spin 
direction (minority-spin) behaves like an insulator with 
Ep inside a band gap. This exotic behavior has inspired 
research not only in the field of basic science, but also for 
applications in spintronics, since the extreme spin polar- 
ization suggests that HMFs are ideal for inducing and 
manipulating the transport of spin-polarized electrons. 

Since HMFs were originally introduced in 1983 0, 
their properties have been explored extensively. Theo- 
retical studies have focused mainly on their ground-state 
properties: the magnetic moments and the origin of the 
gap Q. The stability of the gap was studied with re- 
spect to spin-orbit coupling 3], to surface and inter- 
face states 0, to the presence of defects and to the 
appearance of non-quasiparticle states In parallel, 
the extreme spin polarization has been verified experi- 
mentally in a few compounds Q. However, it is clear 
that the ideal half-metallic property (P = 100%) can- 
not be present at elevated temperatures. Fluctuations of 
magnetic moments will mix the two spin channels, and 
at latest at the Curie point, To the spin polarization 
will vanish together with the magnetization. Thus, for 
application purposes one seeks HMFs with Tq signifi- 
cantly higher than room temperature, reasonably assum- 
ing that the temperature dependence of the spin polariza- 
tion, P(T), approximately follows the magnetization, i.e., 
P(T) cx M(T) 8]. But the theoretical study of P{T) is 
far from trivial. Materials specific, first-principles calcu- 
lations based on density- functional theory (DFT) , which 
capture the physics of hybridizations and bonding es- 
sential to the half- metallic property Q, are designed in 
principle for the ground state and not for excited state 
properties. Nevertheless, adiabatic spin dynamics can be 
approximated within DFT, with successful applications 
in the prediction of To 0> E3] lately also for Heusler 
alloys 0, [H El- 



In many cases, HMFs have more than one magnetic 
atom per unit cell. For example, NiMnSb, a half- Heusler 
compound and prototypical example of all HMF has two 
magnetic atoms: Mn and Ni. It is the d-d hybridization 
between the Mn and Ni minority states that opens the 
half-metallic gap 2|, at least at T = 0. At higher temper- 
atures, directional fluctuations of local moments reduce 
the magnetization. In one-component systems it is well- 
known that the magnetic configuration at each instant 
shows some degree of short-range order: small regions 
present almost collinear magnetic moments, with a local 
spin quantization axis ei oc not necessarily parallel to the 
average moment (the global axis e g i b), while the low- 
energy, long-wavelength fluctuations are more significant 
for the decrease of the magnetization. This behavior con- 
tinues up to and even above To However, rather little 
is known for multi-component systems, when the sublat- 
tices are coupled with different strength so that they can 
lose the magnetic order at different temperatures. 

Motivated by these considerations, we follow two ap- 
proaches to the excited states, both based on DFT (one 
in connection to an extended Heisenberg model, and one 
within mean-field theory) , in order to elucidate the prob- 
lem of magnetization and spin polarization at T > 0. 
Our focus is on NiMnSb. We conclude that the hypothe- 
sis P(T) cx M (T) is not valid. Instead, P{T) can fall off 
much faster than M(T). Furthermore, for NiMnSb we 
find that the Ni moment disorders already at very low 
temperatures, resulting in a susceptibility peak at 50 K. 

We expect that the decrease of P in the presence 
of non-collinear fluctuations arises due to three mecha- 
nisms: (1) Firstly, we have a globally non-collinear effect: 
the local axis e\ oc of a region with short-range order is in 
general not parallel e g i b- Thus, there is always a projec- 
tion of locally spin-up states to the globally spin-down 
direction in the gap region. (2) Secondly, a locally non- 
collinear effect: the short-range order is not perfect, since 
the spin axis of each atom varies with respect to that 
of its neighbors. This effect can be more significant in 
multi-component systems. It is most important if, e.g., a 
spin-up localized d resonance determining the local frame 
is at Ep. The result of this effect is that, even in a lo- 
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FIG. 1: (a): Modification of spin-down DOS around Ef re- 
ducing the spin polarization due to mechanisms (l)-(3) as 
discussed in text, (b): DOS of NiMnSb in the ground state 
(full curve) and in the state when the Ni moment has been 
constrained to an angle of 60° with respect to the Mn mo- 
ment (dashed curve). Inset: Dashed curve: the state when 
Mn moments of neighboring layers along [001] direction form 
an angle of 50° with respect to each other, while the Ni mo- 
ments are allowed to relax. DOS are projected to the global 
frame. 



cal spin frame of reference, spin-up wavefunctions of each 
atom are partly projected into the spin-down states of its 
neighbors within the gap, so that P is diminished. In the 
presence of only the first effect, P(T) is expected to fol- 
low the average magnetization M(T); the second effect 
arises in addition to the first, and P{T) is even further 
reduced. (3) Importantly, there is a hybridization effect 
(change in hybridization strength) , which has been over- 
looked so-far, leading to a closure of the gap by shifting 
the conduction and valence band edges. In Fig. ^ we 
show schematically how the three mechanisms affect the 
minority-spin DOS. Naturally, changes in the majority 
DOS also occur and charge neutrality is conserved. 

We turn now to the specific system, namely NiMnSb. 
We find by first-principles calculations that, in the 
ground state, the Mn local magnetic moment is large 
(Mmii = 3.7 [Ib), while the Ni local moment is much 
smaller (Mm — 0.26 [Ib) and originates mainly from the 
hybridization of the d orbitals of Mn with the ones of 
Ni (the same mechanism which opens the half-metallic 
gap 

0)- 

In particular, the d-d hybridization causes some 
transfer of weight from Mn to Ni in the unoccupied states 
and vice-versa in the occupied ones. Thus, the Ni mo- 
ment is formed, while the Mn moment is reduced. Conse- 
quently, the Ni moment does not really arise from strong 
intra-atomic exchange interactions (Hund's rule), and its 
formation lies at a much lower energy scale than that of 
the Mn moment. We verify this conclusion by perform- 
ing constrained DFT 0] calculations: when the mag- 
netic moment of Mn is constrained by a magnetic field 
to form an angle to its Mn neighbors, it changes very lit- 
tle, even in an artificial antiferromagnetic configuration. 
Contrary to this, the Ni moment vanishes already when 
it is constrained at an angle of 90° with respect to the 



Mn moment, or when the Mn atoms are placed in an 
artificial antiferromagnetic order. 

We approach the finite-temperature properties start- 
ing with a standard recipe: the adiabatic approxima- 
tion for the calculation of magnon spectra |l5|. We em- 
ploy the full-potential linearized augmented plane-wave 
method (FLAPW) as implemented in the FLEUR code 
within the generalized gradient approximation |l7j ] to 
DFT. Total-energy calculations of frozen magnons, per- 
formed at the equilibrium lattice constant of 5.9 A, on 
a dense mesh of 2745 wavevectors q and 4096 /c-points, 
provide the magnon dispersion relations E(q). A sub- 
sequent Fourier transform yields the real-space exchange 
constants Jy between sites i and j (€ {Ni,Mn}) |9|, map- 
ping the system to a classical Heisenberg model. Ther- 
modynamic quantities such as the magnetization curve 
M(T), susceptibility x(T), and Curie point Tc can be 
found within this model by a Monte Carlo method 
In multi-component systems we also consider the sub- 
lattice susceptibility \n and magnetization M n (n is a 
sublattice index). A peak of Xn(T) signals the release of 
the corresponding degrees of freedom, i.e., the sublattice 
magnetization is randomized to a great extent; the total 
susceptibility x(T) presents a peak at To 

While the Mn moment can be treated within the 
Heisenberg model as having a fixed absolute value and 
fluctuating only in its direction, for the weak Ni moment 
the longitudinal fluctuations are energetically as relevant 
as the transversal ones. These considerations require an 
extension of the traditional Heisenberg model allowing 
fluctuations of the magnitude of Mm- This is possible, 
since the energy needed for constraining the Ni moment 
can be calculated within DFT and fitted well by a fourth- 
order function. Within this approximation, at each Ni 
site i the neighboring atoms, placed at sites j, act as an 
exchange field E>i = J2j JijMj polarizing the Ni atom. 
Thus, the energy expression for the Ni atom at i includes 
the magnitude of the local moment Mi and the neighbor- 
induced polarizing field, and the total Hamiltonian reads 



H = - 



iG{Ni} 



+ bM?-B l -M l ). (1) 



The constants a = 18.4 mRy/zxg and b = 42.6 niRy//t|, 
which are fitted to the ab initio total energy results, are 
both positive, giving an energy minimum at Mj = if 
the neighbors have a zero net contribution, i.e., if E>i = 0. 
The second sum of Eq. has to be applied for all Ni 
sites i entering the Monte Carlo calculation, on top of the 
usual first part, which is just the Heisenberg expression. 
Note that, although the Ni moment is small and resulting 
from hybridization, a remnant of intra-atomic exchange 
still exists in Ni, giving an enhanced on-site susceptibility. 
This assists the local moment formation and is reflected 
in the values of a and b. Thus, the system can be regarded 
as an alloy of a strongly magnetic subsystem (Mn) with 



a paramagnetic subsystem with Stoner-enhanced suscep- 
tibility (Ni). 

After calculating the exchange parameters Jy (our re- 
sults agree with those of Ref. [llj), the Monte Carlo cal- 
culation according to Eq. yields the sublattice mag- 
netizations Mmu(T) and Mm(T), and the susceptibili- 
ties XMn{T) and XNi(T), shown in Fig. Evidently, 
the overall thermodynamics are governed by the mighty 
Mn moment. The phase transition is clearly seen by the 
peak in the susceptibility x(T), which grossly coincides 
with the Mn sublattice susceptibility XMn(T). A value 
of Tc ~ 860 K is deduced, verified also by the method 
of cumulant expansion Qs, and lies between the value 
of 940 K calculated in 13] and the experimental value 
of Tc = 730 K. The surprising feature, however, is the 
behavior of Mm(T) and xm(T). Already at low tem- 
peratures, around 50 K, Mm(T) shows a rapid drop and 
XNi(T) a corresponding narrow peak. If we exclude the 
longitudinal fluctuations and work within the traditional 
Heisenberg model, an unpronounced behavior can be seen 
(broad maximum in Xm nsv ( T ) at around T = 300 K 
in Fig. These results show that the Ni sublattice 

magnetic order is lost to its great extent. This behavior 
is traced back to the comparatively weak exchange con- 
stants Jij of the Ni moments to the neighboring atoms, 
as we found by the ab initio calculations. 

It is well-known that the classical Heisenberg model 
cannot capture the very low-energy spectrum of the 
quantum Heisenberg model, therefore the low-T behavior 
of M is usually not well reproduced. However, since the 
longitudinal fluctuations (essential to our model) on the 
Ni sublattice allow for a high rate of energy absorption 
around the crossover temperature (50 K), the classical 
treatment of the Ni sublattice is applicable already at 
such low T. We believe the peak of %Ni and the drop of 
MNi at low T to be connected to the so-far unexplained 
experimental findings of an anomaly in the temperature 
dependence of the magnetization and the resistivity at 
approximately 80 K [IjJ. Assuming that the magnetic 
moments on the Ni sublattice are disordered, the sys- 
tem can absorb energy at a higher rate and this would 
lead to a higher dissipation and subsequent increase in 
resistance. 

The loss of short-range order in the Ni sublattice at 
such low temperatures suggests, at first sight, that half- 
metallicity is easily lost, due to the mechanism (2) of 
local non-collinearity. However, this is not completely 
correct. Since the gap originates from the d-d hybridiza- 
tion between Ni and Mn spin-down states, and since the 
Ni d states are below Ep for both spin directions, a ro- 
tation of the Ni moment causes a hybridization of the 
Mn d spin-down states partly with the Ni spin-up states 
and partly with the Ni spin-down states. Consequently, 
the d-d hybridization remains, but the gap width is re- 
duced. The relevant DOS is shown in Fig. 2Jj, for the 
case of the Mn moments remaining in a ferromagnetic 
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FIG. 2: (color online) (a): Monte Carlo results for the In- 
dependent magnetic properties. The Ni sublattice moment 
Mm (magnified by a factor 10) drops fast at 50 K, fluctu- 
ating transversely and longitudinally; a susceptibility peak 
(XNi(long.+tr.)) is evident. If only transverse fluctuations 
are allowed, XNi(transv.) has a less pronounced maximum, 
(b): Polarization at Ef as function of total spin moment for 
the DLM states for NiMnJ.^MniSb calculated with the CPA. 
The degree x of disorder increases from left (a; = 0%) to right 
(x = 50%), simulating the increase of T from to Tc. 



configuration, while the Ni moments are constrained to 
a 60° angle relative to Mn (and as a result their mag- 
nitude is reduced to 0.16 /is). The reduction of the 
gap width follows from a reduction of the hybridization 
strength: as the Ni moment is reduced, its spin-down 
states move lower in energy and hybridize less with the 
Mn spin-down states. Thus, the spin-down "conduction 
band" ultimately reaches Ep and half-metallicity is lost. 
This is the signature of the hybridization mechanism (3), 
and comes on top of the non-collinear behavior of the 
Mn atoms (mechanisms (1) and (2)). The latter should 
be weak at such low T, since the average Mn moment 
should still be high, dictated by Bloch's T 3/2 law. At 
higher temperatures the locally non-collinear mechanism 
(2) should appear for the Mn moments. In the inset of 
Fig. we show the state when Mn moments of neigh- 
boring layers form an angle of 50° with respect to each 
other. Although mechanism (2) produces a mild effect of 
a finite spin-down DOS at Ep, a protrusion appears just 
0.2 eV under Ep. This is once more related to mech- 
anism (3), triggered by the non-collinear configuration, 
and can cause a collapse of P if the protrusion reaches 
Ep. 

The results so-far do not contain a quantitative esti- 
mate of the polarization as function of T or of the order 
parameter M. Such an estimate requires knowledge of 
the response of the electronic structure to the increase 
of T . We achieve this by proceeding via mean-field the- 
ory, with the disordered local moment (DLM) state 01 
representing the system at T > in a mean-field man- 
ner. A Mn site can then have a "down" orientation 
Mn^ (opposite to the magnetization) with probability 
x, or an "up" orientation Mn^ (parallel to the mag- 
netization). The ferromagnetic ground state (T = 0) 



4 



corresponds to i = 0, while the Curie point corre- 
sponds to x = 0.5. The electronic structure of the 
DLM state for each x is found within the coherent poten- 
tial approximation (CPA), utilizing the Korringa-Kohn- 
Rostoker full-potential Green function method within 
DFT. This yields the average magnetic moment M(x) 
(approximately, M(x) — M(0)(1 — 2x) so x is an alter- 
native order parameter), and the polarization P(x). 

Considering the above, we study the DLM states, 
NiMnj^Mn^Sb, with varying concentration < x < 
0.5. Within this approach, P(T) cannot be directly 
found, if the explicit connection between x and T (or 
M(T)) is not known. However, one can check the hy- 
pothesis P(T) ~ M(T) 8]. As shown in Fig. EJd, such 
a relation does not hold. Instead, from a certain M on, 
P drops much faster than M . The reason can be traced 
back to the change in hybridization as x (or T) increases: 
each Ni atom has on average 4(1 — x) Mn^ neighbors and 
Ax Mn^ neighbors. The occupied d states of the latter 
hybridize with the occupied spin-down states of Ni, push 
them higher in energy and diminish the gap; Ep reaches 
finally the valence band and P collapses. Furthermore, 
an approximation of M(T) consistent with mean-field 
theory can be found by the use of the Brillouin function 
Bj. Doing this for j = 5/2, we assigned the temper- 
ature values shown in Fig. |2t>- Under this assumption 
we see that, up to T = 0.42 Tc (this is about room 
temperature), P{Ep) remains close to 100%, but then it 
drops fast, e.g., P « 35% at T = 0.67 T c . The globally 
non-collinear mechanism (1), not captured by the CPA, 
should be present on top of this behavior. Therefore, the 
initial plateau of P (up to 0.42 Tc) should be corrected 
towards a linear drop as P{T) = M(T)/M(0). 

In summary, we have investigated the behavior of half- 
metallic ferromagnets at elevated temperatures, with em- 
phasis on the properties of the gap and the spin polar- 
ization at Ep. We introduced an extended Heisenberg 
model treating transversal as well as longitudinal mag- 
netic fluctuations, to cover a little investigated situation: 
multicomponent magnets which include subsystems of 
large moments coexisting with paramagnetic subsystems 
exhibiting a Stoner-enhanced susceptibility. The param- 
eters entering the model were determined from first prin- 
ciples. We also estimated the polarization at Ep within a 
CPA averaging. As prototypical system we have chosen 
the half-Heusler compound NiMnSb. Our conclusions 
are the following: (i) In NiMnSb, the Ni sublattice is 
weaker coupled than the Mn one; longitudinal fluctua- 
tions of the Ni moment are energetically as important as 
transversal fluctuations. This leads to an early crossover 
behavior of the magnetization at T w 0.06 Tc where the 
average Ni moment is lost, explaining previous experi- 
ments, (ii) The hybridization creating the gap is still 
present but the fluctuations change its strength. At this 
stage the gap-width is reduced, (iii) At higher T, fluc- 
tuations of the Mn moments introduce a low DOS into 



the gap leading to a mild reduction of P, and, after a 
point, the hybridization changes so much that Ep crosses 
the band edges. This is when the polarization collapses 
(around 0.42 Tc within mean-field theory). The behavior 
of P shown here seems rather general for half- and full- 
Heusler alloys exhibiting half-metallic ferromagnetism as 
we found by additional calculations. E.g., for Co2MnSi, 
the polarization was close to 100% till 0.27 Tc, then it 
droped fast, changed sign at 0.63 Tc and went back to 
zero at Tc. The decisive factor for the thermal collapse 
of polarization is the change in hybridizations due to the 
moment fluctuations. This work calls for experimental 
efforts to measure the sublattice magnetization and the 
spin polarization at Ep as function of T for half-metals. 

We are grateful to P. H. Dederichs for enlightening 
discussions, and to H. Ebert and V. Popescu for providing 
us with their CPA algorithm [2fj| . 
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